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Abstract 

We extended McMillan's Green's function method to study the equilibrium 
spin current (ESC) in a ferromagnet /ferromagnet (FM/FM) tunnelling junc- 
tion, in which the magnetic moments in both FM electrodes are not collinear. 
The single-electron Green's function of the junction system is directly con- 
structed from the elements of the scattering matrix which can be obtained 
by matching wavefunctions at boundaries. The ESC is found to be deter- 
mined only by the Andreev-type reflection amplitudes as in the Josephson 
effect. The obtained expression of ESC is an exact result and at the strong 
barrier limit gives the same explanation for the origin of ESC as the linear 
response theory, that is, ESC comes from the exchange coupling between the 
magnetic moments of the two FM electrodes, J x h,.. In the weak barrier 
region, ESC cannot form spontaneously in a noncollinear FM/FM junction 

when there is no tunneling barrier between the two FM electrodes. 
Pacs numbers: 72.10 -d, 72.25.Dc, 73.23.-b 
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I. INTRODUCTION 



Spin related transport in magnetic hybrid systems has been studied intensively for last 
two decades and considerable progress has been achieved since the discovery of the tun- 
neling magnetoresistance (TMR) effect [1] in ferromagnet (FM) tunnelling junctions. This 
effect originates from the electrical resistances of these junction being dependent strongly 
on whether the moments of adjacent magnetic layers are parallel or antiparallel. The reason 
is that the tunnelling electrons are scattered more strongly in an antiparallel magnetic con- 
figuration than in a parallel configuration. As a result, tunnelling junctions with moments 
in adjacent magnetic layers aligned antiparallelly have larger overall resistance than when 
the adjacent moments are aligned parallel, giving rise to TMR. The reverse effect of TMR 
is the spin transfer effect predicted independently by Slonczewski [2] and Berger [3] about a 
decade ago, in which a sufficiently large spin-polarized current injected from a normal metal 
(NM) into a FM layer can lead to magnetic moment reversal in the FM layer. 

Spin transfer torque occurs in magnetic multilayers with noncoUinear moments and is 
due to the non-conservation of spin current through the interface between NM and FM. 
Due to the presence of noncoUinearity, the component of spin current transverse to the 
magnetization of the layers is not transmitted across the interfaces between NM and FM; 
in other words, the discontinuity of spin current at the interface is the origin of the spin 
transfer effect that can result in magnetization precession or reversal. [4-8] A large number of 
experiments have observed this spin transfer effect in magnetic multilayer structures, [9-11] 
in which one of the FM layers is very thick and works as a polarizer of electric current with 
its fixed moment, and the moment of the thin FM layer may switch when a strong polarized 
electric current perpendicular to the layer plane flows through the layers. 

As the two magnetization directions in the FM/FM tunnelling junction are misaligned, 
an equilibrium spin current (ESC) can flow in the junction without any bias, [7] this phe- 
nomena is analogous to the Josephson effect that the superconductor macroscopic phase dif- 
ference between the two sides of a junction drives a supercurrent through the junction. The 
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dissipationless ESC will disappear if two moments in the adjacent FM layers are coUinear, 
parallel or antiparallel. The existence of ESC J has been verified by many authors using the 
linear response theory, [12-16] and explained as the result of the exchange coupling between 
two magnetic moments and h;, J ~ h/ x h^, Thus the magnetization phase difference 
between the FM electrodes induces ESC in the FM junction as in the Josephson effect. In 
the case of a thin barrier between FMs or the strong coupling limit, the behavior of ESC is 
still unknown in the literature and worth studying. 

In this paper, we study the ESC flowing in a FM/FM junction with two fixed non- 
coUinear magnetic moment using a simple quantum mechanical approach. The present work 
is motivated by two reasons. Firstly, dissipationless spin current in last several years has 
drawn considerable interest and also incurred much debate, such as the controversies in the 
spin Hall effect. [17-20] The study of ESC in FM/FM junctions may shed light on the spin 
current in spin-orbit coupled systems. Secondly, the well-known result of ESC J ~ h; x 
in FM/FM junctions in the hterature was obtained by using the hnear response theory 
[12-16] or heuristic derivation, [21] thus an exact result of ESC is in order. By extending 
McMillan's Green's function method which was originally employed to study the Josephson 
effect in superconductor junction, [22] we obtained an expression of ESC in the FM/FM 
junction at equilibrium. ESC is determined only by the Andreev-type reflection amplitudes, 
[8,23] which is similar to the Josephson effect. [24,25] It was also found that the exact ESC 
expression is slightly different from the linear response result; in other words, the high-order 
tunnelling processes are negligible and the linear response approximation has a wider range 
of application. When there is no barrier between the two FM electrodes, the ESC was found 
to vanish indicating that ESC is a pure quantum mechanical effect. 

This paper is organized in the following way. In Section II, a FM/FM junction model is 
given and the single electron Green's function of the junction system is constructed from the 
scattering coefficients. An analytic expression of ESC in the FM/FM junction is obtained 
in the third section and some discussions are presented. A conclusion is drawn in the last 
section. 

3 



II. GREEN'S FUNCTION 



The FM/FM tunnelling junction is depicted schematically in Fig. 1 and the layer between 
the FM electrodes can be either an insulator barrier or a normal metal. A simple free electron 
model is adopted to describe the FM/FM tunnelling junction and the Hamiltonian reads 

n = — + V{x) - 9{-x)hr(T - 9{x - L)h^-(T, (1) 

where rUe is the effective electron mass and assumed to be identical in all three regions, 
the two FM electrodes and the middle region, and hj. are the internal molecular fields 
of the left and right FM, respectively, cr denotes the Pauli spin operator, and 6{x) is the 
step function. The potential energy V{x) may take different values in different regions but 
remains constant in both FMs. Here the molecular fields and (magnetization with 
energy unit) are not coUinear in our consideration and assumed to be fixed by an external 
magnetic field or other methods. Without loss of generality, we take the spin quantization 
axis of the system to be parallel to the magnetization of the left FM and the direction of 
of the right FM is described by the polar coordinate {9, 0). 

In a free electron model, the energy dispersions of the two FMs can readily be solved 
and they are spin-dependent. In the left FM, the eigenvalue is E± — |^ ± /i; + [// and the 
spinor is the eigenfunction of u^, and in the right FM, E± = |^ ± hr + Ur and the spin 



eigenfunctions are ip'. 
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are the different potential energy in the left (right) FM, k is the wavevector of electron and 



where Ui{Ur) 



± is the spin index. The spatial eigenfunctions are plane waves. Due to spin splitting, there 
are four incoming wave functions with their corresponding outgoing wave functions in the 
left and right FM, as schematically shown in Fig. 2. ^i(x) and ^2(2^) are the wavefunctions 
of the minority spin and majority spin electron injecting from the left FM, rcspctively, while 
^3(0;) and '^i{x) are those injecting from the right FM. For example, the wavefunction of 
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the first type of scattering event ^i{x) is given by 
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,x > L 



In this equation, k''^^^ — ^2m{E — Ui{r) ^ hi{r))/h^ — A;^^ is the spin-dependent wavevector 
along the x-direction in the left (right) FM electrode, E is the single-electron energy which 
is conserved when electrons tunnel through the junction, k// is the wavevector parallel 
to the interface between the different regions and assumed to be conserved in the quantum 
tunnelling process. Nevertheless, this is not a required condition for the following derivation, 
^'i describes a minority-spin electron coming from the left FM lead and being scattered in the 
middle region; the scattering coefficients ai, 61, Ci and di correspond to the normal refiection, 
Andreev-type refiection, [8] transmission without branch crossing, and transmission with 
branch crossing, [26] respectively. For the other 3 scattering processes as well as their 
coefficients a^, fcj, q, and di as shown in Fig. 2 have the same meaning . It is noted that we 
have omitted the parallel plane wave component e*'*//''^^'^^ in the above equation. When we 
use the continuity of the derivatives of wavefucntions at the interfaces to determine these 
coefficients ( Oj, 6^, q, and o?i), the parallel momentum k// should be explicitly taken into 
account. 

With four elementary scattering wavefunctions like those in Eq. 2 above, the single- 
electron Green's function of the junction system can be worked out by the McMillan for- 
mula, [22] which has been further developed by Kashiwaya and Tanaka. [25] This formula has 
been used to treat the Josephoson current in a superconductor junction and relate directly 
the super current to the Andreev reflection amplitudes. The Green's function G^{x,x') 
is proportional to the direct product of the left-going wavefunctions (processes i = 3, 4 
in Fig. 2) and the right-going wavefunctions (processes i = 1, 2), C' ~ ^ l{x)^ r{x') for 
X < x' and ~ r{x) * ^l(x') for x > x', where the hat 'A' denotes the conjugate 



process to the elementary scattering one shown in Fig. 2, and these conjugate scattering 
wavefunctions can be obtained by determining the Hermitian conjugate of only the spinor 
part of the wavefunction not including the spatial part of the wavefunction. The reflection 
and transmission coefficients di, bi, Ci, and di in four conjugate processes have the relations 
di{(f)) ^ ai{-(f)), bi{(f)) = bi{-(f)), Ci{(f)) = and = di{-(f)) (i = 1 ■ ■ - 4), where is 

the azimuthal angle of the magnetization of the right FM. With these scattering wavefunc- 
tions of the elementary processes as well as their conjugate processes, the Green's function 
is then constructed in a linear combination as 



G'{x,x',E) 



ai'^3{x)i/i{x') + a2'if3{x)'^2{x') + Q;3*4(a:)*i(x') + a4'if4{x)if2{x'),x < x' 

(3) 

Here G^{x,x',E) is implicitly a function of the parallel momentum k//. The prefactors ctj 
and Pi {i — 1 ■ ■ - A) can be determined by the boundary conditions that the Green's function 
fulfills, 



e 



^0^ 



(4) 



OX ox % 

G''{x,x\E%=^,+Q = G''{x,x\E%=x,^Q. (5) 

With these two equations, we can directly solve the prefactors ctj and (5i in Eq. 3, which are 
independent of the spatial position x. After some direct algebra, they read 



C3C4 — ^3^4 ' C3C4 — (i3(i4 ' 

-z^ds Z-C3 



C3C4 — 0404 C3C4 — 0404 

Pi = — P2 



C3C4 - dzdi - d2,di 

-z+dz Z-C3 , . 

P3 = — P4 = — J--, (6) 

C3C4 — 0(30(4 C3C4 — 0(3084 

where z± = j^p-- In these solutions, we have employed some detailed balance conditions to 
facilitate our derivation, such as aj(0) = aj(— 0) {i — 1- ■ - A) and /cL&i = k'\b2- Substituting 
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these prefactors into Eq. 3, we obtained the Green's function in the left FM electrode 

{x,x' < 0) as 



G'ix^x'^E) 



y G21 G22 j 



Gil — Z-v exp(i/i;^,.|a; — x'\) + aiz^ exp(— i/c^,.x — ik^_^x'), 

G12 — Z-b2 exp{—ik''_^_x — ik''_x'), 
G21 — z+bi exp{—ik''_^x' — ik''_x), 

G22 — Z- eKp{ik'^_\x — x'\) + a2Z^ exp{—ik''_x — iU_x'). 



(7) 



The Green's functions in other regions such as the middle region between the two FMs can 
also be constructed in a similar manner according to Eq. 3 with solutions given in Eq. 6. 
Although the Green's function is obtained here for the FM/FM junction, this method can 
also be applied to nonmagnetic junctions, and we believe the derivation procedure is universal 
in mesoscopic junction systems. When the retarded Green's function of the studied system 
is worked out, we can in principle calculate the physical observables that we want. Since we 
focus on an equilibrium tunnelling junction, the lesser Green's function is easily obtained 
by the formula 

G<(x, x\ E) = [G%x, x', E) - G'\x, x\ E)] f{E) (8) 

where x' , E) is the advanced Green's function and f{E) is the Fermi-Dirac distribution 
function. The lesser Green's function is related to the spectral function of the electron and 
is very useful for calculating some physical quantities as shown later. 

III. ESC IN THE FM/FM JUNCTION 

In this section, wc first present a general formula of the spin current in an FM and then 
utilize the Green's function obtained in Sec. II to work out the ESC in a FM/FM tunnelling 
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junction. Then as an example, we calculated ESC of a simple FM/FM junction with a delta 
function insulator barrier between the two FM electrodes. 

The local spin density s(r, t) at position r and time t is defined as 

s{r,t) = <if\r,t)kl{r,t), (9) 

where ^(r, t) is a two-component wavefunction and s = |a with a denoting the Pauli spin 
matrices. Taking the time-derivative of s(r, t) and using the Schrodinger equation for a 
single FM, the continuity equation of spin current density is given by 

— + V-J, + S = 0, (10) 

where the respective spin current density and the source term S are 

J, = - — r^t^-v* - (v*)V*i , 

4me^ 



S = ^t(o-xh)*, (11) 

where h is the magnetization of the FM leads (for example, hi, h^.). is the usual definition 
of spin current density and S is referred to as the spin source term or spin torque that leads 
to spin rotation when the spin is not coUinear with the magnetization h. Therefore this 
spin source term gives rise to a spin current flow which is not included in J^. Using the 
lesser Green's function in Eq. 8, G^{xt,x't') — i{^\x't')^{xt)) where (• • •) is the quantum 
statistical average, the spin current density and the source term S above in a steady state 
can be rewritten as 



Arue 



S = lim / i^Tr [(h X a)G<{x', x)\ , (13) 

where the trace is taken over the spin space. Substituting the Green's function given by 
Eq. 7 and Eq. 8 into the above equation, we obtain the ESC expression in the FM/FM 
junction as 
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'^s^-J^Y.^'^ {(^i - + 1:^ ^M-Kk'+ + /(E) (14b) 

k// 

= 0. (14c) 

Here the summation over A;// includes the contribution of all transverse modes to the spin 
current, where /c// is assumed to be conserved in the quantum tunnelling process. According 
to the above equation, the spin current density is determined by the Andreev-type reflection 
coefficients hi and 62; this result is very similar to the formula of the Josephson current which 
is directly related to the Andreev reflection coefficients. [24,25] The spin current density 
in Eq. 14 is exponentially x-dependent so that the integral over energy and transverse 
density of states will make it disappear in bulk FM not far away from interface. This 
phenomena is same as the spin transfer effect that the transmitted spin current is also 
dependent on position and will vanish after summing all possible transverse modes with 
different wavevectors. Stiles and Zagwill [6] estimated the characteristic length of this spatial 
spin precession was about l/Zc/ {kf Fermi wavevector) because only the electrons near the 
Fermi energy contribute to the electric current. For our case in Eq. 14, all the band electrons 
contribute to ESC so that Js is expected to decay much more quickly with distance from 
interface. 

Apart from the spin current density J^, which can be calculated directly from its def- 
inition, additional spin current flow arises from the source term S according to Eq. 10. 
According to the Gaussian theorem, another spin current J^oc can be deflned from this 
source term as 

4:l^j's^dx = I ^5^i?e|(A;?,-^L)(^-^)(l-expH(fc^ + fcL)x])}/(£;) (15a) 

k// 

Jl.^ jy^dx^-j ^ J] 7rn|(A;^-fcL)(^ + ^)(l-expH(fc^ + fcL)x])} /(£;). (15b) 

kj/ 



The 2;-component of Jsoc is zero. The x-dependent part of Jgoc is just the — in Eq. 14, so 
that the sum of Jgoc and is independent of x and keep constant from the interface to the 
bulk FM, which does not disappear away from the interface, 

k// 

When a spin current enters the bulk FM, it will be absorbed by the lattice in FM because 
the magnetization h rotates the spin and the sum of J^oc and therefore remains constant. 

As an example, we calculated a simple case of the FM/FM tunnelling junction, the 
potential energy in both FM is Ui — Ur — 0, and the insulator barrier is described by a 
delta barrier UoS{x) where Uq denotes the strength of the barrier. The magnitudes of the 
magnetizations in two FM are equal hi = hr- Using simple quantum mechanics method, we 
obtain the Andreev-type reflection amplitudes as 

61 = 2{k_ — k+)k+ sin 9 exp{i(p)/ A, 
62 = 2(A;_ — k+)k^ sin 9 exp(—i(f))/ A, 
A = k'i + kl + 6k+k- - 2kl + AikQ{k+ + kJ) - {k+ - kJf cos 9, (17) 

where k± = k^^"* and /cq = mUo/fi^. Substituting these quantities into Eq. 16, the total ESC 
is then given by 



J 



/ ^EMlM(^)](^+-M/(^)h/xh, (18) 



where h;(0, 0, 1) and hr(sin 9 cos 0, sin 9 sin (f), cos 9) are just the unit vectors of the magnetiza- 
tions in the left and right FM. The ESC depends on the directions of the two magnetizations 
not only through the cross product term h; x but also on the quantity A in Eq. 17. In 
the finite barrier case, we can neglect the cos^ term in A so that Eq. 18 can reproduce 
the result obtained with the linear response theory, i.e., the exchange coupling between two 
magnetization leading to ESC. At the opposite limit, zero barrier between two FMs Uq — 0, 
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Eq. 18 indicates there is no ESC flowing in the junction. This is the pure quantum mechan- 
ics effect, the reflected spin direction by a barrier rotates with respect to the incident spin 
direction and an imaginary part can appear in the quantity A of Eq. 17, which was regarded 
as [6] one of reasons that spin current is not conserved at the interface when a polarized 
charge current flows through a NM/FM junction. For ESC discussed here, the reflections 
by a barrier leading to an additional phase to the incident wavefunction which is necessary 
for the formation of ESC in the noncoUinear FM/FM junction. It is quite different from the 
Josephson effect between two superconductors, in which a cooper pair carries supercurrent 
through the junction even with zero barrier. 



IV. SUMMARY 

We have studied ESC flowing in the noncoUinear FM/FM tunneUing junction by extend- 
ing McMillan Green's function method. The single-electron Green's function of the junction 
can be constructed by the scattering coefficients, reflection and transmission amplitudes 
which can be directly calculated by simple quantum mechanics method. In the derived 
formula of spin current density, we found that the exact result of ESC does not have a big 
difference from that obtained by the linear response theory, in other words, the high-ordered 
tunneUing process contribute httle to ESC. It was also found that when there is no barrier 
between two FMs, the ESC will disappear for the reflected spin state from a finite barrier 
has an additional phase which is crucial for the formation of ESC. 
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FIGURES 

FIG. 1. Schematic of a FM/FM tunnelling junction with two noncoUinear magnetizatons in 
the left and right FM electrodes. The middle region M at < x < L can be either an insulator 
barrier or a normal metal. The z-direction is taken as the spin quantization axis parallel to the 
left magnetization h;. 

FIG. 2. A schematic picture of four elementary scattering events. Every injecting electron 
will be scattered into four outgoing wave functions, two reflections and two transmissions, e.g., 
^'i(a;) denotes a minority spin electron k\. in the left FM injecting into the middle region with two 
reflection (ai and bi) and two transmission (ci and di) wavef unctions. 
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Figure 1 wang and chan 




Figure 2 wang and chan 
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